MPI-Ph/93-77 
TUM-T31-50/93 

November 1993 



Theoretical Uncertainties and Phenomenological 
Aspects of B ^ Xsj Decay* 

A. J. Buras/'2) M. Misiak,i)t M. Miinzi) and S. Pokorski2)t 

Physik- Department 
Technische Universitdt Miinchen 
D-85747 Garching, Germany 

Max- Planck- Institut fiir Physik 
Werner-Heisenberg-Institut 
Fdhringer Ring 6 
D-80805 Miinchen, Germany 



ABSTRACT 

We analyze uncertainties in the theoretical prediction for the inclusive branching ratio 
BR[B Xg'j]. We find that the dominant uncertainty in the leading order expression comes 
from its /i-dependence. We discuss a next-to-leading order calculation of -B ^ Xgj in general 
terms and check to what extent the /^-dependence can be reduced in such a calculation. 
We present constraints on the Standard and Two-Higgs-Doublet Model parameters coming 
from the measurement of 6 — 57 decay. The current theoretical uncertainties do not allow 
one to definitively restrict the Standard Model parameters much beyond the limits coming 
from other experiments. The bounds on the Two-Higgs-Doublet Model remain very strong, 
though significantly weaker than the ones present in the recent literature. In the Two-Higgs- 
Doublet Model case, the b 57, Z ^ bb and b — > crUr processes are enough to give the 
most restrictive bounds in the Mh± — tanj3 plane. 
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1. Introduction 

The b ^ S'y decay is known to be extremely sensitive to the structure of fundamental 
interactions at the electroweak scale. As any Flavour Changing Neutral Current (FCNC) 
process, it does not arise at the tree level in the Standard Model (SM). The one-loop W- 
exchange diagrams that generate this decay at the lowest order in the SM (fig. |I|) are small 
enough to be comparable to possible nonstandard contributions (charged scalar exchanges, 
SUSY one loop diagrams, Wr exchanges in the L-R symmetric models, etc.). 
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Figure 1: Leading order one loop diagrams generating the 6 — S7 transition in the SM. 



Among all FCNC processes, the 6 — >■ S7 decay is particularly interesting because its rate 
is of order G'j^aQED, while most of the other FCNC processes involving leptons or photons are 
of order Gl^a^^^. The long-range strong interactions are expected to play a minor role in the 
inclusive B Xg'y decay.0 This is because the mass of the b-quark is much larger than the 
QCD scale A. Moreover, the only relevant intermediate hadronic states ip^s are expected 
to give very small contributions, as long as we assume no interference between short- and 
long-distance terms in the inclusive rate Therefore, it has become quite common to use 
the following approximate equality to estimate the B —>■ Xg'j rate 



T[B ^ XcCUe] T[b ^ cei>e 

where the quantities on the r.h.s are calculated in the spectator model corrected for short- 
distance QCD effects. The normalization to the semileptonic rate is usually introduced in 
order to cancel the uncertainties due to the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements and factors of m\ in the r.h.s. of eq. (|l]). 

As indicated above, this ratio depends only on mt and in the Standard Model. In the 
extensions of the SM, additional parameters are present. They have been commonly denoted 
by ^. The main point to be stressed here is that R is a calculable function of its parameters 
in the framework of a renormalization group improved perturbation theory. Consequently, 
the decay in question is particularly suited for the tests of the SM and its extensions. 

One of the main difficulties in analyzing the inclusive B — > decay is calculating 
the short-distance QCD effects due to hard gluon exchanges between the quark lines of the 
leading one- loop electroweak diagrams. These effects are known [0-0 to enhance the 6 — > 37 
rate in the SM by 2-5 times, depending on the top quark mass. So the 6 ^ 57 decay appears 
to be the only known process in the SM that is dominated by two-loop contributions.^ 

With the above conclusion in mind, one can realize how astonishing it is that the 6 — > 57 
decay is already measured! The recent CLEO report [|1^] gives the following branching ratio 
for the exclusive B — > K*-^ decay 



BR[B K*j] = (4.5 ± 1.5 ± 0.9) x 10"^ (2) 

^Xs can be identified by requiring that E-y > {ra\ — Tnjj)/2mB 
■^apart from otlier q — > q'7 transitions which are very hard to measure 
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Therefore, we can say that the inclusive B Xg'j branching ratio is measured to be larger 
than the one in eq. @, and smaller than 



BR[B X,7] < 5.4 X 10"^ {at 95% C.L.). (3) 

These experimental findings are in the ball park of the SM expectations based on the 
leading logarithmic approximation. It is then not surprising that after announcement of the 
results (H) and (H) several analyses appeared in the literature with the purpose to restrict 
the parameters of the SM or put bounds on its extensions. 

In view of these new developments, it is the purpose of the present paper to reanalyze 
the theoretical uncertainties in the calculation of the ratio We find, in agreement with 
ref. [0], that the dominant uncertainty in the existing leading logarithmic calculations is 
due to choice of the renormalization scale /i. Such an uncertainty, inherent in any finite 
order of perturbation theory, has been recently analyzed in several FCNC processes, such as 
B^-B^ and K^-K^ mixing or rare K- and B-decays |13|, |T^. It has been demonstrated 



in these papers, that the inclusion of next-to-leading order corrections reduces considerably 
the /i-dependence of the relevant amplitudes. 

Since B — >• is dominated by QCD effects, it is not surprising that the scale-uncertainty 
in the leading order calculation of this decay rate is particularly large — it amounts to around 
±25%. It follows that the restrictions on the SM or its extensions which can be obtained 
with help of the experimental findings (§) and @) and the leading order approximation are 
substantially weaker than found by other authors. 

Unfortunately, it will take some time before the //-dependences present m B ^^7 can 



be reduced in the same manner as it was done for the other decays []T2|, |T^, |T^. As we 
will describe in the following, a full next-to-leading calculation of 5 — > Xg'^ would require 
calculation of three-loop mixings between certain effective operators. Before one undertakes 
such an effort, it is instructive to make a formal analysis of the considered decay at the 
next-to-leading level and to check to what extent the /i-dependence can be reduced once all 
the necessary calculations have been performed. 

Our paper is organized as follows: In section 2 we summarize the results of the leading 
logarithmic calculations. Subsequently, we discuss various uncertainties present in the exist- 
ing formulae. In section 3 we analyze the restrictions from eq. on the parameters of the 
Standard and Two-Higgs-Doublet Models, taking into account the uncertainties discussed in 
the previous section. In section 4 we generalize the renormalization group analysis beyond 
the leading logarithmic approximation, we list the calculations which have to be done (or 
have been already done), and demonstrate how the //-dependence will be reduced. Section 
5 gives a brief summary of the paper. 

2. Present status and theoretical uncertainties 

Let us start by summarizing the results of the leading logarithmic calculations 
In a compact form, they can be written as follows: 

^^I|^_^6^ (4) 

where 
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with z = = as(Miy)/as(/i), and 



g{z) = l-8z^ + 8/ - z^ - 2AzHn z. (6) 

The function g{z) is the phase space factor in the semileptonic b-decay.0 The numbers hi 
and Oj are given and discussed in the appendix. Here we only mention that the sum of all 
the hi vanishes. Finally, the Standard Model values of the coefficients Cf'\My/) are Wb\ 



Cf\Mw) = l (7) 



A{x - ly 24(x - 1)3 ^ ' 

^(^^) = i(^3TjI^^^ + 8(0: -1)3 

2 

where x = Some details of the derivation of eqs. (|)-(||), as well as the reason for 

introducing the symbol Cf^'^^^jj) will be given in section 4. 

The renormalization scale ^ present in eqs. and has to be of order m^, but need 
not to be exactly equal to nii,. The QCD coupling constant ^^(yu) at any renormalization 
scale /i can be expressed (as it has recently become popular) in terms of its value at = Mz- 
The leading logarithmic expression is just ag^jj) = as{Mz)/v{fi), where 

vW = l-f3o — ^ — ln\—j (10) 

and Po = 11 — |/. In our case, the number of active flavors / equals 5. 

Before starting the discussion of uncertainties, let us illustrate the relative numerical 
importance of the three terms in the expression (^ for cf^^^^ . For instance, for irit = 
130GeV, 11 = 5GeV and a,(Mz) = 0.12 one obtains 



C)"^^^^(/i) = 0.689 Q"^(Mvk) + 0.087 C^iMw) - 0.161 C^iMw) = 

= 0.689 (-0.161) + 0.087 (-0.086) - 0.161 =-0.280. (11) 

In the absence of QCD we would have Cf'^'^^^^i) = C^^\Mw) (in that case one has rj = 1). 
Therefore, the dominant term in the above expression (the one proportional to C2°'*(Mvy)) is 
the additive QCD correction that causes the enormous QCD-enhancement of the b ^ s'j rate. 
It originates solely from the two-loop diagrams. On the other hand, the multiplicative QCD 
correction (the factor 0.689 above) tends to suppress the rate, but fails in the competition 
with the additive contributions. 

The equations (|l]) and (^-(^ summarize the complete leading logarithmic (LL) expres- 
sion for the B —>■ Xsj rate in the SM. Their important property is that they are exactly 
the same in many interesting extensions of the SM, such as the Two-Higgs-Doublet Model 
(2IIDM) or the Minimal Supersymmetric Standard Model (MSSM).| The only quantities 
that change are the coefficients ci^\Mw), Ci^\Mw) and C^^^Mw) (see refs. [|, for 
details). 

^Note, that at this stage we do not include the 0{as) corrections to r(6 cev) since they are part of 
the next-to-leading effects discussed in section 4. 

^However, in a general SU{2)l x SU{i)ji x U{1) model they get modified, because some new "effective 
operators" enter - see ref. il6| . 
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Let us now have a critical look at eqs. (|l]) and (|D-(0)) ^i^^ theoretical and exper- 
imental uncertainties present in the prediction for BR[i? —>■ Xs'y ] that can be made with 
help of these equations. 



(i) First of all, eq. (|I]) is based on the spectator model. For a long time there has been a 
question of how good this model is and what is its relation to QCD. Fortunately, during 
the past two years, this relation has been better understood. The spectator model has 
been shown to correspond to the leading order approximation in 1/mb expansion. The next 
corrections appear at the 0{l/ml) level. The latter terms have been studied by several 
authors [jl8| with the result that they affect BR[i? Xg'-f ] and BR[i? XcCUe] by only 
a few percent, as long as one uses = m^°''^ in the leading term. In view of much larger 
uncertainties discussed below, the exact size of the (9(l/m^) corrections is not essential. The 
effects due to the difference between m^"'"^ and mb{fi = mf,) are next-to-leading QCD effects 
which we discuss in point (v) below, and in section 4. In our phenomenological analysis 
presented in the next section, we will assume the error due to the use of spectator model 
in B ^ Xg'-f to amount to ±10%. This error is also meant to include interference between 
the short-distance contributions and contributions from the ipXg intermediate states [|l|. It 
could be significantly larger, had interference terms in all the exclusive channels the same 
signs, but we can see no reason for such a correlation. 



(ii) As already mentioned, the normalization of the B Xg'y rate to the semileptonic 
rate in eq. (|I|) has been introduced in order to cancel the factors of mf on the r.h.s. of 
this equation. It also drastically reduces uncertainties due to the CKM matrix elements (see 
point (iii) below). The price to pay is an additional uncertainty coming from the ratio z = ^ 
in the phase-space factor for the semileptonic decay. However, this ratio is experimentally 
known better than the actual masses of both quarks. And this is the only combination in 
which these two masses enter our expressions (P and (1^)-([10|). 

The phase-space factor g{z) changes rather fast with z around -2 = |- In the next section 
we will use z = 0.316 ± 0.013 |T^. The resulting error in the ratio R is then around 6%. 
The actual accuracy of determining z can of course be subject to discussion. However, even 
if one assumes that the error due to z is twice larger than 6%, it still remains significantly 
smaller than the ones we discuss below. 



(iii) The error due to the ratio of the CKM parameters in eq. (^ is small. Assuming unitarity 
of the 3 X 3 CKM matrix and imposing the constraint from the CP-violating parameter 6x 
we find ^ 

^^^^^ = 0.95 ± 0.04 (12) 

The quoted error corresponds to 0.036 < |Kfe| < 0.047 0], \Vub/Vcb\ = 0.08 ± 0.02 and 
Bx = 0.7 ± 0.2. Eq. (|12|) has been obtained for nit set to 200 GeV. For smaller values 
of rrit the central value is practically the same, and the error decreases down to 0.03 for 
TTit = 100 GeV. 



(iv) There exists an uncertainty due to the determination of a^. This uncertainty is not 
small because of the importance of QCD corrections in the considered decay. Now all the 
extractions of which give conservatively the range 0.11 < as{Mz) < 0.13 in the 
MS scheme, include next-to-leading order corrections both in the expressions for and 
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in the relevant processes used for as determination. Unfortunately, the next-to-leading 
order corrections to B ^ Xg'j are unknown at present. So for this decay it is not possible 
to consistently use as{Mz), say in MS scheme. For this reason, the usage of two- loop 
expressions for t] in eq. @ can certainly be questioned. However, if one takes the same 
initial condition as{Mz) = aj^{Mz) for the evolution of «s(/^), then the difference between 
i?'s in eq. (^) resulting from using the two-loop and the one-loop expressions for ag amounts 
only to 4-9%. On the other hand, the difference between the ratios R of eq. (|) obtained with 
help of aj^{Mz) = 0.11 and 0.13, respectively, is between 15 and 22%. So the inaccuracy 
in determination of a<j(M^) is more important here than the NLL effects in the evolution 
of a<;(/i). In the next section we will use just the leading logarithmic expressions for and 
include only errors due to varying as{Mz) between 0.11 and 0.13. 

(v) The dominant uncertainty in eq. (P comes from the unknown next-to-leading order 
contributions. This uncertainty is best signaled by the strong //-dependence of the leading 
order expression (^), which is shown by the solid line in fig. 0, for the case rrit = 150 GeV. 
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Figure 2: /i-dependence of the theoretical prediction for the ratio R, for mt=150 GeV and 
as{Mz) = 0.12. The solid line corresponds to the leading order prediction. The dashed lines 
describe possible next-to-leading results (see section 4). 

One can see that when /x is varied by a factor of 2 in both directions around ruf, ~ 5 GeV, 
the ratio (|) changes by around ±25%, i.e. the ratios R obtained for /i = 2.5 GeV and 
H = 10 GeV differ by a factor of 1.6. This is in a rough agreement with what has been found 
In section 4 we will argue why varying /i by a factor of 2 in both directions 



in ref. [11 



properly estimates the size of unknown next-to-leading QCD corrections. 

The dashed lines in fig. ^ show the expected /^-dependence of the ratio (Q) once a complete 
next-to-leading calculation is performed. The /t-dependence is then much weaker, but until 
one performs the calculation explicitly one cannot say which of the dashed curves is the 
proper one. The way the dashed lines are obtained is described in section 4. 

For the purpose of the next section, we will estimate the size of the unknown next-to- 
leading terms by varying /i by a factor of 2 in both directions around ^ 5 GeV in the 
existing leading logarithmic expressions (i)-®- 

There is one more comment we should make at this point: The scale used in the numerator 
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of the ratio rj = as{Mw) / ctsifJ') (the "matching" scale) does not have to be exactly equal 
Mw but rather can also be varied by a factor of 2 around M^. The resulting changes in 
the ratio R are then roughly twice smaller than the ones we have found by varying the low- 
energy scale around mi,. However, this uncertainty is also due to next-to-leading effects. So 
we will treat it as already taken into account in the estimate we have made by varying the 
low-energy scale, i.e. we will not vary the matching scale in our analysis described in the 
next section. 

(vi) Finally, we have to mention that there exists a 4.6% error in determining BR[B Xs'~f] 
from eq. (|l]), which is due to the error in the experimental measurement of BR[B XcCUe] = 
(10.7 ±0.5)% p. 

We treat the mj-dependence in eqs. (|D-(^ as an "interesting uncertainty", i.e. all our 
predictions will be given as functions of rrit. We will never add this uncertainty to the other 
errors. 

The above discussion of uncertainties is to a large extent model- independent. None of the 
errors listed in points (i)-(vi) was due to the coefficients ci°\Mw), Cj°^(Mvi/) and ^(M^/) 
(except for NLL contributions to them).0 Therefore, in all extensions of the SM where only 
cI°\Mw) get modified (hke in the 2HDM or the MSSM), the above estimates remain valid. 
Of course, in extensions of the SM, some additional uncertainties in cf°''(Mw) may occur. 
Apart from the dependence of cI^^\Mw) on the parameters of these models (like e.g. tanjS 
or charged scalar mass in the 2HDM), which is a welcome feature, we can also encounter 
additional uncertainties due to the SM parameters: For instance, Cf°''(Miy) in the MSSM 
contain ratios of the ordinary and SUSY-sector mixing angles (see ref. [0). 

3. Phenomenological consequences for the Standard and Two-Higgs-Doublet 
Models 

In the present section we describe in what manner the inclusion of theoretical uncer- 
tainties affects the limits on the Standard and Two-Higgs-Doublet Model parameters that 
can be obtained from B —>■ Xgj measurement. 

Fig. § presents the SM prediction for the inclusive B Xs'~f branching ratio. In order 
to obtain the error bars presented in this figure, we have added in quadratures the errors 
(i)-(vi) listed in the previous section. Adding in quadratures does not mean, of course, 
that we can attribute any statistical meaning to these errors. However, it is the best one 
can do at present in order to estimate the total uncertainty in the theoretical prediction for 
BR[i? — i> ]. By adding in quadratures we want to take into account that there is no 
correlation between the errors (i)-(vi). 

Comparing the theoretical SM prediction with the experimental results in fig. |^, one 
can conclude that it is rather improbable that the future measurements will contradict 
this prediction, in view of its large uncertainty. Comparison with experiment can become 
sharper only by making the theoretical prediction more precise — first of all by calculating 
the complete NLL short-distance effects. 

One can also ask a related but somewhat different question of whether, assuming the 
present theoretical uncertainty, a future precise measurement of the inclusive B Xg'y rate 

^Thc effect of the QCD evolution between the rrit and Mw scales can be treated as NLL contribution as 
long as ln{mt/AIw) < 1, i.e. when rrit < 220GeV. This contribution has been already calculated in ref.psf 
and found to give around 10% effect. 
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Figure 3: Predictions for B — >■ Xgj in the SM, as a function of the top quark-mass, with the 
theoretical uncertainties taken into account (see the text). 



may tell us something about the SM parameters. 

Fig. ^ presents the lower limits on the top-quark mass following from the assumption 
that the future experimental result will satisfy the inequality 

BRtheo^B ^ X 7I + 2 X f theoretical \ ^ ( measured \ 

^ * I error bar j \ decay rate j 

for a certain value of m^. Such mt is then treated as disallowed. The three curves correspond 

to different possible results of the measurement of BR[i? — > ] (3, 4, and 5 times 10""^). 
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Figure 4: Lower limits for rrit in the SM as a function of as{Mz) (see the text). 

If we trust our estimation of the theoretical uncertainty, the conclusion is that we may 
get a useful lower bound on the top quark mass provided the branching ratio is close to the 
present experimental upper bound. It has to be stressed, however, that the limits in fig. ^ 
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are quite sensitive to the treatment of the errors. Even if one only symmetrized our errors 
in fig. H, these hmits would move up by around 10%. 

One could also ask what upper limits on the top quark mass (as functions of as{Mz)) 
can be obtained by assuming BR[i? — ^ Xg'j ] is measured below the theoretical prediction, 
and the inequality 

BR^'-[B X,7] - 2 X f theoretical \ ^ f experimental \ 

\ error bar / \ upper bound / 

is satisfied. We find that upper bounds for rrit that lie below 200 GeV can be obtained 
only if one assumes that the inclusive rate will be measured smaller than around 1.7 x 10~^. 
This is not very encouraging, in view of the present measurement quoted in eq. (0), and the 
expectation that the B K*'y rate forms 10-20% of the inclusive rate [0, depending 
on the model used for making such an estimate. 

Let us now discuss the limits on the Two-Higgs-Doublet Model. 



The relevant part of the 2HDM lagrangian |^ which is responsible for the nonstandard 



contributions to 6 — > 57 has the following form: 



AG 



F 



XULVcKAiMnDn + YURMuVcKMDL + h.c. (15) 



Here is the charged scalar field, U and D represent the column vectors of the up- and 
down-quarks, respectively, Vckm is the CKM matrix, and Mu {Md) is the diagonal mass 
matrix for the up- (down-) quarks. The quantities X and Y are given by tanjS = V2/V1, i.e. 
the ratio of the vacuum expectation values of the two scalar doublets. In the most popular 
version of the 2HDM, usually called "Model 11" , the masses of the down-quarks (up-quarks) 
are proportional to Vi (f2)-[] Then one has 

1 Vo 

X = — = tanp = — . (16) 
Y vi 



Figure 5: Charged scalar exchanges contributing to 6 — > 57 in the 2HDM. 

The charged scalar exchanges shown in fig. ^ create extra contributions to the coefficients 
ci'^\Mw) and cI^\Mw) (see ref. §): 



AC^°^(Mvk) 



3y' -2y -5y^ + 3y 1 
-My + — — + 



Q{y - If ^ 12{y - If tan^(3 



Sy' - 2y' 



■Iny + 



ly^ - 5y^ + 7y 
12{y - 1)3 



(17) 



-y , -y^ + ^y . i 

Iny + — — + 



2(2/ -1 



A{y - 1)2 tan^/3 



-r 



Aiy-lf 



Iny + 



-y^ + by^ + 2y 
24(y-l)3 



''The notation in rcf.j|] is opposite. 
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where y = mf/Mjj±. 

As we have aheady said, equations @-(|D remain unaltered. So from eq. ([I7D we 
can immediately obtain the 2HDM predictions for the B — > branching ratio. They 
are presented in fig. |], for the case tanjS = 2 which is large enough to make the terms 
proportional to 1/ tain? [3 in eq. (|T^ practically irrelevant. For smaller tan/5, the results grow 
rapidly as tan/? decreases. 
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Figure 6: Predictions for B in the 2HDM, as a function of the top-quark mass, with 

the theoretical uncertainties taken into account (see the text). 

For most of the top-quark masses, the results for Mh± = lOOGeV in fig. ^ are in 
contradiction with present measurements, even if the theoretical uncertainties are taken into 
account. On the other hand, comparing figs. ^ and ^ we can conclude that in the Mh± = 
GOOGeV case we are not able to distinguish between the SM and the 2HDM, independently 
of the accuracy of future measurements (but assuming the present accuracy of the theoretical 
predictions). 

Fig. [^presents the limits in the tan(3-MH± plane, obtained by assuming that the excluded 
values of tanf] and Mh± are those for which the inequality (|l^) holds. The presented limits 
are found with help of the experimental input from eq. (^. The three solid curves represent 
lower bounds on Mh± obtained from b ^ S'-f for three different values of the top quark 
mass: 120, 150 and 210 GeV. The dashed curves show bounds coming from the Z — > 66 
and b — * crur decays. The shown Z ^ bb bound is based on the results of ref. and 
corresponds to rrit = 150 GeV. The b — * crUr bound is given by the simple formula 



tan/? < 0.54: 



M 



H 



IGeV ' ' 

The b — >■ ctUt- and Z ^ bb bounds are complementary to the b ^ s'j bound. All other bounds 
in the tanj3-MH± plane (e.g. the ones coming from B^-B^ mixing and other FCNC processes 
considered in refs. [^) are weaker than the combination of the three above-mentioned ones. 

Our b s'-f bounds are much weaker than the ones found in refs. |30, 31], because we 
have taken into account the uncertainties in the theoretical prediction. If we did not take 
them into account, our bound for rrit = 210 GeV would lie as high as the dotted curve in 
fig. 0. 
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As one can see, the b ^ s'-f bounds in fig. |^ are practically independent of tan(3 for its 
values larger than 2. From the theoretical standpoint, large values of tan(3 are interesting 
since tanjS ~ mf/m^ would naturally explain the large splitting between these two masses 



without necessity of introducing order-of- magnitude different Yukawa couplings |3^. For 
these tanjS^s it is convenient to plot the tan/?-independent bounds in the mt-MH± plane. 
They are presented in fig. M. 
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Figure 8: The b —>■ sj lower limits for the charged Higgs boson mass as functions of rrit for 
tan/3 = 2. See the text for the description of the curves. 

The solid line in this figure shows the lower bound for Mh± obtained using the requirement 
described in eq. (p!^). If we require the theoretical prediction for BR[i? Xg'y ] be only 
one, not two error bars away from the experimental upper limit, then our lower bound on 
Mh± is described by the dashed curve in fig. ^j. The dotted line in this figure corresponds to 
assuming that the theoretical prediction is given by its central value, with no uncertainties 
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whatsoever. We present the three curves in order to demonstrate the relevance of taking the 
theoretical uncertainties into account. 



Our limits are much weaker than found in ref. |3^, and they do not exclude the t H^b 
decay of the top-quark. However, they still have to be considered very strong. For the most 
interesting values of tanjS they are still more restrictive than the bounds from the Z ^ bb 



decay ||2^, |3^ , as well as the bounds from other FCNC processes . The expected new data 
on B ^ Xs'j in the nearest future will further improve the bounds on the 2HDM, but the 
possible improvement is limited by the theoretical uncertainties. Comparison between figs. 
§ and 1^ shows that for Mh± larger than around 600 GeV, the present theoretical inaccuracy 
does not allow to distinguish the 2HDM and SM predictions. 

4. The structure of B — > Xs'j beyond leading logarithms 

In this section we describe a complete next-to-leading calculation of i? — > Xgj in general 
terms, and check to what extent the theoretical uncertainties due to the /x-dependence can 
be reduced after performing such a calculation. 

The QCD corrections to the B Xs'j decay contain large logarithms ln{M^/ml) which 
have to be resummed with help of the Renormalization Group Equations (RGE). In order to 
do this, one introduces an effective hamiltonian built out of operators of dimension higher 
than 4. 

Heff = -Vi,l^,:^|:g.(/i)a(/i) = -Vt,V*^Q^{^,)C{f,) (19) 

where Vij are the elements of the CKM matrix, Qi{n) are the relevant operators and Cj(yu) 
are the corresponding Wilson coefficients. The complete set of operators necessary in the 
6 — > S7 case is the following [Q: 



V-A 



Ql = (s„C^)v-a(C/3&' 
Q2 = iSaC")v-AiC(}b^)v-A 

Q3= {sJ)°')v-AY.q{.W^)v-A 

Qa= iSab'^)v~AEq{m°')v-A 



(20) 



V+A 



Qe = {sJ}'^)v--AY.q\qi3T)V+A 

Qs = 8f^mfes„a^-(l + -f5){TT0bf'G% 

where {qq')v±A = Qlfii^ i 75)9' and a^i, = |[7/x,7i/]- Qi,2 are the current-current operators, 
Qs-e the QCD penguin operators, and Qt.s the "magnetic penguin" operators. In the latter 
operators we have neglected the terms proportional to the s-quark mass since they give only 
a correction of order m^/m^ to the decay rate. 

The inclusion of the couplings in the definition of Q7 and Qs allows one to discuss the 
renormalization group evolution in the full analogy with the case of the usual AS = 1 or 
AB = 1 hamiltonians for nonleptonic decays. 

Following p3 we write 



Cifi) = U{fi,Mw)C{Mw) (21) 
where the renormalization group evolution matrix is given generally by 
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f/(m..,„,) = r,e.p/^^^_^^d9'^ (22) 

with mi < Here Tg denotes such ordering in the couphng constants that they increase 
from right to left. Next, 7(5') is the 8x8 anomalous dimension matrix of the operators Qi, 
and I3{g) is the usual renormalization group function which governs the evolution of ag- 
Keeping the first two terms in the expansions for '^{g) and (3{g) 



2 4 
\ _ ^ (0) 9 _|_ ^(1) 9 

^^^>~^ 167r2^ (167r2)2 



167r2 



2\2 ' 



(23) 



where (3q = II - f / and pi = 102 - f /, one finds Q: 



U{mi,m2) 



1 + 



An 



4tt 



J 



(24) 



where tj^^\mi,m2) denotes the evolution matrix in the leading logarithmic approximation 
and J summarizes the next-to-leading corrections to this evolution. If 



(25) 



where 7^5''' denotes a diagonal matrix whose diagonal elements are the components of the 
vector 7'-°^ then 

Ui'\mi,m2) = v\ (H^]\ with 
For the matrix J one gets 



^a^(mi) / 



7 



(0) 



D 



2Po 



(26) 



where the elements of S are given by 



J = vsv~ 



(0) A 



G 



Wo 2Po + 



7: 



(0) .,(0) 



(27) 



(2^ 



7 



with 7!°'' denoting the elements of 7^°-' and Gij the elements of G. 



The ratio as{m2) / as{mi) in eq. ( pBD has now to be calculated with help of the two-loop 
renormalization group equation for as(/i). Once as{Mz) is treated as an initial condition for 
this equation, the corresponding solution takes the form 



1 - 



(3i asjMz) In v{fi) 
Po 4:71 v{n) 



(29) 



where v(fi) is exactly as in eq. (|T0|). 

To the same level of accuracy we expand C{Mw) as follows 



w 



(30) 



The nonvanishing Cf\M\Y) have been given in eqs. (0)-(|^) 
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The leading logarithmic approximation to Ci{fi) is completely known. It is obtained by 
setting J = and C^^^ = in eqs. (|2^) and (|30|) respectively. In the appendix, we explicitly 
give the matrix 7^°^ calculated in a certain regularization scheme called the HV scheme. 
The scheme dependence of 7*^°-' will be discussed below. Here we would like only to remark 
that the 6x6 submatrix of 7'-°^ describing mixing of (Qi, Q2, Qe) and the 2x2 submatrix 
for {Qj, Qs) follow from one-loop calculations. On the other hand, the leading order mixing 
between these two sectors results from two-loop diagrams. 



Let us enumerate what has been already calculated in the literature and which calculations 
are still required in order to complete the next-to-leading calculation of Cj(/i). Beyond the 
leading order, 7^^^ and C^-^\Mw) have to be calculated. Actually, the 6x6 submatrix of 
7^^-' describing the two-loop mixing of (Qi, ■■■,Qe) and the corresponding initial conditions 
m C^^^Mw) have been already calculated in refs. |3^, The remaining ingredients of a 
next-to-leading analysis of C(/i) are: 

(i) The two-loop mixing in the {Q7,Qs) sector of 7'-^^ This calculation is currently be- 
ing performed ||36| . 



(ii) The three-loop mixing between the sectors {Qi, ■■■yQe) and (QtjQs) which with our 
normalizations contribute to t"*-^-*. 



(iii) The 0{as) corrections to Ct^M^) and Cs{M]y) in eqs. (||) and (|). This requires eval- 
uation of two-loop penguin diagrams with internal W and top quark masses and a proper 
matching with the effective five-quark theory. An attempt to calculate the necessary two- 
loop Standard Model diagrams has been recently made in ref. The effective theory 
one-loop diagrams with Q7 and Qs insertions have been calculated in ref. p4[. However, the 
finite parts of the effective theory two-loop diagrams with the insertions of the four-quark 
operators (see figs. 5 and 6 of ref. 0) are still unknown. 



All these calculations are very involved, and the necessary three-loop calculation is a 
truly formidable task! Yet, as will be evident from the following discussion, all these calcula- 
tions have to be done if we want to reduce the theoretical uncertainties in 6 — > 57 to around 
10%. 

It is important to stress that a next-to-leading calculation performed without resumming 
large logarithms ln{M^/ml) would not be more accurate than the present leading-order one 
where these logarithms are resummed. This is why the calculation of the three-loop mixing 
is unavoidable. 



Once the evolution of the coefficients Ct(/i) down to the scale /i ~ is performed, 
one has to evaluate the matrix elements of all the operators at this scale. 

In the leading order, only the tree-level matrix element of Qj and the one-loop matrix 
elements of Qi — Qq that are of order have to be included. The latter matrix elements 
vanish for the on-shell photon in any 4-dimensional regularization scheme and also in the HV 
scheme, i.e. dimensional regularization with the 't Hooft-Veltman definition of 75 All 



but the last calculations of the leading order QCD corrections to 6 57 used the NDR 
scheme, i.e. the dimensional regularization with fully anticommuting 75 .0 In this scheme, 
the one-loop 6 — *• 57 matrix elements of some of the four-quark operators Qi, ...,Qg do not 

® A complete calculation using the "dimensional reduction" scheme is not yet existent - see ref. p9] . 
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vanish for the on-shell photon but are proportional to the tree-level matrix element of Qi 
(even for the quarks off-shell). 

The regularization scheme-dependence of leading order matrix elements is a peculiar 
feature of the 6 ^ 57 decay. It can arise because one- loop mixing between the (Qi, Qe) 
and (Q7, Qs) sectors vanishes. In consequence, what usually would be a next-to-leading 
order effect is only a leading one. And what would usually be next-next-to-leading (like the 
above-mentioned three-loop mixing) is only next-to-leading. 

This peculiarity causes, that at the leading order it is convenient to introduce the so- 
called "effective coefficients" for the operators Qj and Q^: First, we observe that in any 
regularization scheme one can write the one-loop matrix elements of the four-quark operators 
Qi) Qe as 

{Qi)one loop — ViiQl) trzey ^ — 1, 6 (31) 

for the on-shell photon but quarks possibly off-shell. In consequence, the leading order 
b ^ s'j matrix element of -ffe// is equal to the tree level matrix element of 

-V,,V,:^Cr"if^)Or (32) 

where 

Cr'^(/.) = Cf)(^) + Ey.Cr )(/.). (33) 

i=l 

In the HV scheme all the ?/j's vanish, while in the NDR scheme y = (0,0,0,0, — |, — 1). 
This regularization scheme dependence is cancelled by a corresponding regularization scheme 
dependence in 7'^'') (see ref. p for details). Consequently, the quantity cf^^^^^) defined in 
eq. ( ^3|) and used in eq. (H) is scheme- independent. 

The numbers Ui in eq. ( |3TD come from divergent, i.e. purely short- distance parts of the 
one-loop integrals. So no reference to the spectator-model is necessary here. This is why 
one is allowed to treat the expression ( |32D as the proper effective hamiltonian that has to be 
inserted in between hadronic states in a calculation of, say, the B K*'y decay. 

An "effective coefficient" for the operator introduced in a similar way equals to 



i=l 

where the numbers Zi are defined by the one-loop b —>■ s gluon on-shell matrix elements of 



; d 

the four-quark operators: 



{Qi)one loop — Zi{Q8)treey ^ — 1,...,6. (35) 

All the Zj's vanish in the HV scheme, while in the NDR scheme we have z = (0, 0, 0, 0, 1, 0). 
For later convenience we introduce a scheme-independent vector 

C(°)^//(/.) = [C?\f,),....,CP{f,),Cr'\f^),Cr'\f^)] ■ (36) 

From the RGE for C^^\fi) it is straighforward to derive the RGE for C^^^^^^ {^i). It has the 
form 

^^j^r''^\^^ = lT''cf^'^\^^ (37) 

where 
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Ifi' + ELi Vklfk - VjlT? - Zjl87 when i = 7 and j = 1, 6 
^T" = \ 7]? + ELi ^J7]? - ^,7^? when . = 8 and j = 1, 6 (38) 
[ "yjf otherwise. 

The matrix scheme-independent quantity. We give it exphcitly in the appendix. 

At the next-to-leading level the situation is not as simple. One has then to include one- 
loop 6 — > S7 matrix elements of Q7 and Qs, and two-loop matrix elements of Qi, Qq. As 
long as the photon and the external quarks are put on-shell, we can still formally write the 
considered matrix elements as proportional to the tree-level matrix element of Q7. One could 
naively expect that the proportionality constants depend only on the ratio yu/m^ because after 
neglection of the s-quark mass, we have only a single scale in the process. Unfortunately, they 
also depend on the infrared regulator used in the calculation.^ However, their /x-dependence 
is still a purely short-distance effect and can be found with help of the existing leading-order 
results. 



As discussed in ref. the /^-dependence of {Q{fi)) can be calculated in a renormaliza- 
tion group improved perturbation theory, as long as fi is not too small. Taking the scale mb 
as a reference point, the /i-dependence of (Q(/u)) is simply given by 

{Q^if^)) = {Q^{mb))U{mb,fi) (39) 



with U given already in eq. (^4]). The scale mi, is chosen here only for convenience. Inserting 
eqs. (^) and (|39D into the final expression for the b —>■ S'j amphtude: 

A{b ^ S7) = -VtbV*^ E(Q.(/x))a(/x) (40) 

V / i=i 

one finds that the /i-dependence present in C(/t) is cancelled by the corresponding /t- 
dependence in {Q{fi)). 

Let us now write the sum of the tree and one-loop on-shell matrix elements of Qr and Qs 
calculated at /t = rrif, as 



{Qi{mb)) 



as{mb) 

^■n H Ti-j 

An 



{Q7{mb))tree, ^ = 7,8, (41) 



and the sum of one- and two- loop on-shell matrix elements of Qi, ...,Qq as 



{Qiinib)) 



asirub) 



{Q7imb))tree, i=l,...,6. (42) 



As long as the infrared divergencies are regulated dimensionally, the r^y's do not contain 
dimensionful parameters, i.e. they are either pure numbers or depend only on the infrared 
regulator e^. 

If /t ~ TUb, there are no large logarithms present in {Q{fi)), and the resummation of 
logarithms is not needed in eq. (|39D. Therefore expanding eq. (|39| ) around /t = m^, we 
can reproduce the /i-dependent terms that would appear in the straightforward perturbative 
calculation of {Q{fi)). Using eqs. (^)-(|4^) and expanding the matrix U{mb,fi) from eq. 



^The infrared divergencies are removed later by including soft gluon brcmsstrahlung - see ref. p4 
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(|39| ) around fi = rrih we find, that the b ^ s'y amphtude calculated in terms of the coefficients 
renormalized at any scale /i ~ rrih can be written in the following form: 

A{b ^ S7) = -VtbV,:^ D {Qj{m,))tree (43) 

where 

D = C',"(,) + {cr'"(i^^'"lr^ + C;°V)r.v) • (44) 

All the Cj's in the expression (|33|) for C^^^ {^) are now taken up to next-to-leading accuracy.|^ 
The term proportional to In"^ comes solely from expanding {Q{ix)) around /i = m;,. The 
tree-level matrix element of Q-j is denoted as {Q7{mi,))tree in order to point out that the mass 
mi, in the normalization of this operator is renormalized at the scale fi = rrih. 
A next-to-leading generalization of eq. (^) takes the following form: 

R = = l^^Fl^r (45) 



(46) 



with D defined in eq. (^4|), and the factor F given by 

/ \ 2 



The quantity k,{z) is a sizable next-to-leading QCD correction to the semileptonic decay 

(47) 



K,{z) ~ 1 — 



which tends to cancel with the other factor in eq. (^ 



/ mfc(/i = mb) \ _ S asjmb) , . 

V mf'^ j 3 TT ' ^ ^ 

and, in consequence, the expression (^) differs from unity by only around 8%. 

The origin of the latter factor is very simple: We have normalized the 6 — >^ S7 rate to the 
semileptonic rate in order to cancel factors ml in front of the spectator-model expressions. 
However, in the semileptonic case, all five factors of mi, come from the on-shell external 
lines, while in the case of 6 — > 37 two of the m^'s come from the normalization of Q-/, and 
should be taken in the same renormalization scheme as the coefficients. This becomes even 
more obvious when one realizes that the anomalous dimension matrix element 777'' = y is 
dominated by the anomalous dimension 7^) = y = 8 of the mass m^ which normalizes Q7. 
If we did not include the running mb{fi) in the operator normalization, then our coefficients 
would run differently than described in eq. (^. 

Let us now return to the complete next-to-leading expressions (^)-(^) for the 6 — *• S7 
amplitude. From eq. (||) and from the RGE (^) for C(°)"^^(/i) one can easily see that the 



dominant /i-dependence in the leading logarithmic coefficient Cf^^^^^) will be cancelled 
by the logarithms coming from the next-to-leading matrix elements. Consequently, the 



This refers to the first term in D. In the term proportional to one can stih keep the leading logarithmic 
(^jo)e//^ Note, that r^y's are multiplied by the usual, not the effective coefficients. The resulting scheme- 
dependences should cancel with the ones present in ri-j and C^^-'(/x). 
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resulting /i-dependence of the obtained result will be one order higher than the accuracy of 
the performed calculation, as it always must happen in any perturbative calculation. 

We can explicitly check the cancellation of /x-dependence, because the coefficients at the 
logarithms in eq. ( ^Ij) are given by the matrix ■y^^^^ff which is already known (see the 
appendix). However, any change of fi is equivalent to a shift in the unknown constant terms 
Tij. Consequently, a meaningful analysis of the /^-dependence must also include rjy. These 
constants are generally renormalization-scheme dependent, and this dependence can only be 
cancelled by calculating 7*^^-' in the same renormalization scheme. Since this point has been 



extensively discussed in several papers ||3J, we will not repeat this discussion here. 



However, it is clear from these remarks, that in order to address the /i-dependence and the 
renormalization-scheme dependence as well as their cancellations, it is necessary to perform 
a complete next-to-leading order analysis of along the lines presented above, and to 
calculate the constants r^r in eqs. (^l]) and (|4^) . 

The program of calculating the complete next-to-leading contributions and reducing the 



/^-dependence has been already accomplished in several other weak decays. In ref. |T2[ the 
/f-dependence in top-quark contributions to K^-K^ and B^-B^ mixings have been considered. 
A similar analysis of the charm quark contribution to K^-K^ mixing has been recently done 
in ref. [0]. Next, in ref. the /t-dependences in rare decays Kl tt'^uu, tt'^uu, 



Kl jjLjl^ B ^ II and B — >■ Xgyv have been studied. In all these papers it has been 
demonstrated explicitly that the inclusion of next-to-leading corrections considerably reduced 
the //.-dependence present in the leading order expressions. The remaining /i-dependence can 
be further reduced only by including still higher order terms in the renormalization group 
improved perturbation theory. 

In view of the tremendous complexity of the next-to-leading calculations for 6 — >■ S7, we 
are not yet in a position to make a similar analysis for this decay. We can however perform 
a simple exercise with the terms we already know, in order to get a rough idea by how much 
a future next-to-leading order calculation could reduce the /i-dependence shown by the solid 
line in fig. ^. 

We proceed as follows: We set J = and C^^^ = in eqs. ( plD and (|30|) respectively, i.e. 



we calculate C'^f^{jj) in the leading logarithmic approximation. Cf'^'^^n) has been already 
given in eq. (13). For the remaining coefficients we find 
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)(/i) = ir^-i - ^r^i+O.OSlOr^o-^ose _ 0.140377^0.4230 _ o.OllSr^-o-s^s^ + 0.005 V'^^^^ 
Cf )(/i) =-|r/-i - j^r/i+0.098V-4086 ^ 0.1214r7-0-4230 + o.0156r/-0-8994 ^ Q^QQ2Qr^o.u5(i 
Cf (/i) = -0.0397r7™ + OMUr]-^-^^^^ - 0.0025r7-o-8994 + 0.0304r/0-i456 

cPiix) = 0.0335r/0-4086 + o.0239r7-o-4230 _ 0M62ri-'^-^^^^ - 0.01 12770- ^^se 

and 

cr-(,).(cf(M.).H||)„ii_ 

- 0.9135r/0-^o«'^ + 0.0873r7'0-^23o _ q^q^yIt]'^-'"''^ - 0.0209r]'^ '^''' 
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Next, we set|^ rn = in eq. (^31), except for for which we consider values —2, and 2. It 
is clear from eqs. (|4l|)-(^2D that all the rj7's can take values of order 1. We choose to vary 

because it is the operator Q2 that has the largest coefficient — the coefficient {rrih) 
is around 1.1. The "next-to- leading" results for the ratio R found this way are presented as 
functions of /i by the dashed lines in fig. ||. 

One can see that the /i-dependence is much weaker than for the leading-order results (the 
solid line), but not completely absent. The initial ±25% uncertainty is reduced to around 
±3%. The dominant /^-dependence caused by factors of order as{'mi,)ln{mb/ fi) cancels out, 
but we still have noncancelling higher order terms like a1{mb)ln'^{mb/ fJ^) or al{mb)ln{mb/ fi). 
We have organized our exercise in such a way that both kinds of the latter terms do appear, 
in spite of that we have set a lot of the next-to-leading corrections to zero. The higher order 
/i-dependence could be, in principle, removed by including still higher order QCD corrections. 

Now it is easy to see that, indeed, a reasonable variation in possible non-logarithmic 
next-to-leading terms, namely changing from —2 to 2, causes a similar change in the final 
results as varying /t by a factor of 2 in the leading order terms. This supports the method we 
have used for estimating the size of the NLL corrections in our phenomenological discussion 
in section 3. 

The simple exercise presented above explicitly illustrates the cancellation of the /i-depen- 
dence present in the leading-order expressions. It cannot, however, replace the full next-to- 
leading order calculation. Indeed, we still do not know which of the dashed curves in fig. |] 
is closest to the true next-to- leading result. Our exercise is, nevertheless, rather encouraging 
as it shows that the theoretical prediction for B —>■ Xgj will be considerably improved once 
all the next-to-leading calculations are completed. 



5. Conclusions 



i) Among many theoretical uncertainties present in the calculation of the B Xg'y 
decay, the most important one is due to the absence of a complete calculation of next-to- 
leading short-distance QCD corrections to this decay. It causes around ±25% inaccuracy in 
the prediction for this decay rate, while all the other uncertainties are around or below 10%. 

ii) Once the present experimental data are taken into account, it is quite improbable 
that a future precise measurement of the B Xg'j rate will contradict the Standard Model 
prediction, in view of its large uncertainty. However, if the future measurement of the 
inclusive rate is close to the present upper bound, then some useful correlation between nit 
and as{Mz) will emerge. 

iii) In spite of the theoretical uncertainties, the present experimental data still put quite 
strong bounds on the parameters of the Two-Higgs-Doublet Model. The b —>■ 57, Z —>■ bb 
and b —>■ crur processes are enough to give the most stringent bounds in the MH±-tanf3 
plane. However, our b s'-f bounds are weaker than found by other authors. 

iv) The theoretical prediction can be significantly improved by making a systematic next- 
to-leading calculation within the framework of renormalization group improved perturbation 
theory. Such a calculation is very involved, since it requires a calculation of finite parts of 
many two-loop diagrams, and divergent parts of three-loop diagrams. However, it has to be 
done if we want to make the uncertainty coming from short-distance effects smaller than the 
other uncertainties. 

v) The uncertainty in the leading order prediction is best signaled by its strong dependence 

Here we assume we have already removed the infrared divergencies from aU the ri^'s by adding soft 
gluon bremsstrahlung. 
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on the renormalization scale /i at which the RGE evolution is terminated. We have described 
in detail how this dependence is reduced in a full next-to-leading result, and performed a 
simple exercise to check it numerically. As a result, we can conclude that once a next-to- 
leading result is found, its /x-dependence should be below 5%, for fi changing from mb/2 to 
2mb. 
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6. Appendix 

Here we give the scheme-independent matrix ■y^^^^J'f defined in eq. (^S|) that governs 
the leading-logarithmic QCD corrections to 6 — > 57 . In the HV scheme this matrix is equal 
to the matrix 7^°'' of eq. (p3l). 
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(0)e// 



-2 


6 

















3 


6 


-2 


2 


2 


2 


2 


416 


70 


9 


3 


9 


3 


81 


27 








22 
9 


22 


4 
9 


4 
3 


464 
81 


W + 3/ 








6-1/ 


-2 + 1/ 


--f 


-f 

3'' 




6 + i/ 














2 


-6 


32 
9 


-T-3/ 








--f 


-f 


--f 


-16 + 1/ 


— -g]-M + -^d 


-4-W/ 




















32 
3 























32 
9 


28 
3 



(49) 



where f = u + d, and u (d) is the number of up- (down-) quarks in the effective theory. In 
our case u = 2 and d = 3. 

The above matrix has been obtained under the assumption that the disagreements be- 
tween refs. 0, and |Q concerning 757'*, 75g'' and 7g8'' in the NDR scheme are resolved in 
favour of ref. 
of ref. 
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. More information about these disagreements is contained in the appendix 
Here we only mention that these disagreements make the numerical leading order 
results differ by only around 1%, so they are irrelevant for our phenomenological discussion. 
Note, that our normalization of Qr and Qs differs from the one used in ref. 0. 

The scheme-independent numbers and hi in eq. (^ are given by the eigenvalues and 
eigenvectors of ■y^^^^^-f .|^ Their explicit values are given by 



= ( if, if, 4, -g, 0.4086, -0.4230, -0.8994, 0.1456) 
hi = i fiff, -|f|, -f, -j^, -0.6494, -0.0380, -0.0186, -0.0057) 

With help of Mathematica all the above numbers can be found exactly. However, 
some of them are so complicated that it is much better to keep them numerically. 



^^according to eqs. (|25|) and (G6h with 7'^"^ replaced by ^^'^^^ff 
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